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A quintessential Inflation (QI) scenario from Mimetic Dark Matter (MDM) is pre¬ 
sented in this paper. This scenario, which is based on the MDM model presented 
by Chamseddine and Mukhanov [1], uses a potential that is defined on three time 
intervalsdhe first during inflation, the second after the end of inflation and 

the third after that. The resulting energy density of the universe is constant during 
inflation, followed by that of a matter/radiation dominated universe, and Anally 
ends with a constant energy density corresponding to dark energy. The scale factor 
has an accelerating expansion nature during and after inflation. It will be shown 
how this is still a viable scenario, even if the scale factor after inflation is not that 
of a decelerating De Sitter universe. 


Chapter 1 
Introduction 


Quintessential Inflation models have been introduced drastically in the literature in 
an attempt to link inflation to the later stages of the universe’s evolution |1]-[8]. The 
key element in this unihcation is the fact that both inflaton, the held describing 
inhation, and quintessence, the held describing dark energy, are both dynamical 
scalar helds that are describing an accelerating expansion of the Universe. However, 
one cannot avoid to mention that some models have used gauge helds to describe 
these phenomena [9], nevertheless we are not going to tackle on these in this paper. 
The way this is usually done is by adding the Lagrangian of a scalar held to that of 
general relativity, while assuming an FRW metric: = dt'^ — a{t)‘^6ijdx^dx^ where 

a{t) is the scale factor and 6ij is the Kronecker delta. This will result in the following 
action: 

S = + !/(.#.)) (1.1) 

where g is the determinant of the metric, R is the Ricci scalar, 0 is the QI held, the 
dot denotes derivative with respect to time and V(0) is the potential which describes 
the dynamics of the held. Instead of introducing scalar helds from outside into the 
Lagrangian, Chamseddine and Mukhanov wrote the physical metric in the following 
wayd]: 

d/ll. = d/ii.(d"^5a0^/30) = Pg^,u (1-2) 


where is an auxiliary metric, 0 is (for the moment) a random scalar held and da 
denotes partial derivative with respect to a;“. In this way, one might say that the 
conformal mode of the metric has been isolated, for the physical metric is invariant 
under a conformal transformation of the auxiliary metric. Moreover, one can see 
from 


1.2 that such an isolation results in the following constraint equation: 




(1.3) 


which will be essential in specifying the scalar held later on. From here, the resulting 
action would be: 
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where is the matter content of the system. Thus, the gravitational field acquires 
an extra longitudinal degree of freedom, in addition to the two transverse ones 
representing the graviton. This extra degree of freedom will prove to be useful. 
Now, the resulting equation of motion (upon varying the action with respect to the 
metric) gives the Einstein tensor in terms of the stress-energy tensor and 
an extra term 

QflU _ ( 15 ) 


where 


t""' = (G -T)g‘“g''^d„ 


( 1 . 6 ) 


with G and T being the trace of the Einstein tensor and the stress-energy tensor, 
respectively. Now compare 1.6 to that of a perfect fluid, = {p + p)u'^u'^ — 
where p is the energy density, p is the pressure and is the 4-velocity which 
satisfies the normalization = 1. By doing the following identihcations: p = 0, 
p = G — T and = g^°‘da4>, one can see that this extra degree of freedom imitates 
“dust”. Therefore, one is not obliged to introduce any type of matter from outside 
to explain the phenomena attributed to dark matter, rather now one has to extract 
hidden fields from the metric to explain these phenomena. Choosing the FRW 


metric, one can see that by solving 1.3 , the resulting scalar field would be: 


= t 


(1.7) 


which will be used through out this paper. 

Furthermore, such isolation can be introduced into the action in the following way 

0 - 0 : 

S = f - 1) - V(4,) + L„| (1.8) 

where A is a Lagrange multiplier. Now, varying the action with respect to the metric, 
we get the “modihed” Einstein equation: 


- 2\d„4id„4i - g^„V = 0 


(1.9) 


Taking the trace of 1.9 , we get the following relation for A: 


A = -[G-T-AV) 


which means, after plugging back in 1.9 


G,, = {G-T- AV)d^<t)d,<t) -I- g^,V{(t)) + T, 


/il/ 


( 1 . 10 ) 


( 1 . 11 ) 


As was done above, the hrst two terms on the right hand side of 1.11| can be 
identified with the stress energy tensor of a perfect fluid with pressure: 


p = -V 


( 1 . 12 ) 


and energy density. 


p = G-T-3V 


(1.13) 
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The time-time component of 1.11 gives the Friedman equation; 


= -~e = — / aVda 


(1.14) 


Multiplying |1.14 by a^, and differentiating with respect to time, while substituting 
y = we get: 


y- = 0 


(1.15) 


(see |3] for detailed derivation). 

From here, the work of this paper is based on choosing an appropriate potential for 
QI, plug it in 1.15 and study the resulting expansion of the Universe and get its 
energy density. We will end up with a conclusion and future works. 
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Chapter 2 
The Model 


In this chapter, we will consider a quintessential inflation model for MDM. The 
structure of the potential was inspired from a paper by Peebles and Vilenkin [S], 
although the exact details are not the same. Moreover, the content of the potential 
is based on the dynamics of a slow rolling held, but using MDM. We will then go 
onto considering the appropriate potential in MDM that would produce almost the 
same effect of QI. 


2.1 Inspiration from Slow Rolling Cosmology 

Consider a slow-rolling scalar held, with the following potential: 

V = 0-“'^ (2.1) 

with 0 = ln(f) |1] is the QI held and a is a constant; thus we are using a power law 
potential. The scale factor that comes from a slowly rolling held (using non-modifled 
General Relativity)is : 


a = aoexp{- 


a 


c-“). 


( 2 . 2 ) 


'3(2-a) 

with oo being a constant of integration. The energy density of such a held would 
be: 


1 


P = 


a 




(2,3) 


where Mpi is the Planck mass. With an appropriate choice of a, this model shows 
an exponential expansion of the universe, but with an energy density that goes like 
This energy density is that of radiation and matter [TU]. On the other hand, if 


-2 


t 

we take the following potential: 
the scale factor then is: 


V = /3e-^ 


a = aoexp(-l^t) 


with the same definition for 0 = ln(f). Moreover, the energy density is now: 

1 


P = 


3M2 




(2.4) 

(2.5) 

( 2 . 6 ) 
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which is a constant. From here, we see that to produce an energy density that repre¬ 
sents matter-radiation dominated universe (i.e goes like t~‘^) directly after inflation, 
and reaches an asymptote, the potential must be a combination between the two. 
Combining the two potentials together, while substituting the form of 0, we get a 
potential of the form: 

V = Ar^ + Bt (2.7) 


So, let’s try to see what physics will be produced from MDM if we use a polynomial 
potential, with 1.7 as the scalar field. Moreover, Peebles and Vilenkin used the 


following potential: 


A(c/)^ + M% 0 < 0 


V = 


AM8 

04 + ’ 


( 2 . 8 ) 


> 0 


where M is a parameter to be fixed by data. Although there are some differences 
between this work and the one in [8] , it would be interesting to see what could such a 
potential (that is one which is defined on two intervals of 0, or equivalently two time 
intervals) produce when used in MDM, and it will be shown that some important 
results appear. 


2.2 The Potential 


Let’s use the following potential in MDM: 


2n/ 1 

V = —(1 - a)(t - ta)^“ + -|a(« - (o)““ + 


(2,9) 


The origin was shifted to tg, the time at which inflation ends. As compared to 2.7 


this has been done in order to make sure that the dynamics are centered at the end 
of inflation . Moreover, the choice of the coefficients is made in such a way that no 


clustering of constants occurs later on. Plugging this potential in 1.7 , we get: 
«(! — a) 


y 


-{t-to) "--[a(t-to)^ “-7 0]^]2/ = O 


the solution of this equation will give us the scale factor to be[TT]: 
a = aoexp[- 


^ \2-a 1 


‘3(2 - aA ~ 


and an energy density for the mimetic matter: 


a 


P = 


(t-ta)' 




( 2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 


One can see that if a is very small, p oc t ^ at the beginning, that is near the end 
of inflation, and then as t —)■ oo, p —>■ 0^. So far, what we have is exactly 
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the behavior we expect. What remains is hxing the parameters a and (3 to produce 
the desired measurable quantities. However, there’s still something wrong with this 
potential. First, the energy density and the scale factor might diverge, unless we 
have a good choice of the parameter a at the boundaries. Second, if t < to, and 
we have a fractional power in the energy density and the scale factor, we will get 
imaginary numbers. This is something definitely we don’t want in real measurable 
quantities. 

Therefore the solution will be as follows: we will separate the potential into two 
parts, one before inflation {t < to) and the other after inflation {t > to). We will 
then match these two values ai t = to- This way, we will have the term to — t during 
inflation {t < to) and the term {t — to) after inflation. By doing this, we have solved 
the issue of having imaginary numbers. Now, concerning the divergence issue, we 
look at the form of the scale factor and the energy density in 2.10 and 2.11 . To 
avoid divergences, during inflation, at t = to, 2 — a must be positive, so must be 

1 — a. Therefore the solution to avoid divergence at f = to as we approach it from 
the left, is to have: 

a < 1 

Now, for t > to, keeping the same form of the potential, our concern is at oo, since 
there we don’t want the energy density to diverge, rather we want it to be a very 
small number. Moreover, the scale factor should not diverge at t = to- Therefore, 

2 — a' > 0 and a' — 1 > 0(we are using a' just to distinguish it from the constant 
during inflation). So for the post-inflation phase: 

1 < a' < 2 


The final result for the potential that would produce a quintessential inflation model 
in MDM is: 


"^(1 ““ ^o) ^ +-[e{t — to) ^ —/3']^, 


t > t 


0 


V = 


- !)((„ - 1)'-^ + 1[(2 - £)(«„ - 1)'-^ + eif, t < t„ 


(2.13) 


and the corresponding scale factor is : 


- to)" - - to)]. 


t > to 


a = < 


aoexp| gp‘_ (to -tf ■ + ^/3(«o - t)], t < h 


(2.14) 


while the energy density: 


P = 


1 . 2-e 


mli\t-to)^ 
1 . 


/S'] 


•n2 


(to-t)' 


ri + 


t > to 


t < to 


(2.15) 
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where e is an infinitesimal nnmber. Now, to determine /3, we have to nse the nnmber 
of e-folds of inflation. If inflation is to last for 70 e-folds, then: 


r-io 


N = 


Hdt = 70 


(2.16) 


where ti is the time at which inflation is snpposed to have started. According to 
the model first presented by Gnth, inflation shonld start at ti = 10“^® and end at 
= 10-32 


t 


Flagging in these nnmbers into 2.15 

X 10^2 


we get: 


(2.17) 


On the other hand, /?' is determined by matching the valne of the energy density at 
infinity to that of the cosmological constant [13]. This will resnlt in 


/?' ^ 7(3) X 10 


-23 


(2.18) 


Before continning into checking the resnlts of the model, there’s one last issne that 
needs to be tackled. It is apparent from the form of the energy density in 2.15 
that it diverges at t = to- This might mean that there’s a discontinnity in the 
energy density at the end of inflation. We can approximately solve this issne by 
looking at how mnch time it takes p to go from cx) to the valne at t = to if we 
are approaching it from the left (i.e. nsing the expression of the energy density for 
t < to). If we ping in the valne of /? in p for t < to, we get the energy density 
at t = to to be of the order of 10^®®. Setting this valne to be that of the field 
for t > to, and calcnlating the time interval, it turns out that it takes the energy 
density approximately 10-®®s to go from oo to 10^®®. Since it is a very short period 
of time, we can insert yet a third interval, directly after inflation, which extends for 
only 10-®®s. During this time, the energy density, the scale factor and the potential 
take approximately constant values corresponding to those at the end of inflation. 
Therefore the third part (now is the second)will end 10-®®s after the end of inflation 
and there will be no divergence in the form of p when taken in t > to + 10“®®, 
and the continuity problem is therefore solved. This maneuver wont affect any 
previous or later calculations, since it extends over an extremely short period of 
time. These equations will result in the plots below for the scale factor and the 
energy density (we have used an e = 0.01). From the first plot, it is clear that the 
scale factor is increasing with a > 0, which implies it is an accelerated expansion 
of the universe. The expansion during inflation is much steeper than that after it, 
which is exactly what’s needed, since the universe cannot keep on accelerating at the 
same rate as that during inflation. Furthermore, concerning the energy density plot, 
the graph shows a constant energy density during inflation, which is a characteristic 
of inflation. In addition, the energy density reaches an asymptote as t ^ oo, which 
is nothing but Quintessence. The inflation parameters, with the above expressions 
of the scale factor, are consistent with the conditions for inflation: 


H 

e = -—=oc 10 


-35 


« 1 ; 


V = 


dine 

~dN~ 


0 


(2.19) 
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Figure 2.1: Plot of the logarithm of the scale factor as a function of the logarithm 
of time, for the two regimes: during inflation (red) and after inflation(blue). It is 
clear that there’s a huge expansion in the universe during inflation, and a moderate 
one afterwards. The value of e that has been used is 0.01 


Figure 2.2: Plot of the logarithm of the energy density as a function of the logarithm 
of time for the two regimes 
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Now we will discuss what type of perturbations does this model lead to. As was 
pointed out in [3j, in order to have a difference between the short wavelength pertur¬ 
bations and the long ones, a term of the form (7 is a constant)must 


be added to the action. Moreover, the short wavelength perturbations have been 
shown to be independent of the choice of the potential. On the other hand, for long 
wavelength perturbations, we do have a dependence on the choice of the potential, 
for it depends on the scale factor. The integral in the fluctuation is taken over the 
period of inflation, that is from ti = 10 “^® to to = 10“^^. At the end of inflation, 
the two terms in exp of 2.14 die away. Therefore we can say that the integral is 
dominated by the lower limit, and since we have an exponential expansion, we can 
approximate the form of the scale factor to be a ~ exp(/9(t — to))- This will make 
the integral much easier to calculate. From here, we get the perturbations in the 
scalar held to be: 


6(f) = A- 
a 


a^dr] 


A 1 
_ _ 

/6 ~ H 


( 2 . 20 ) 


which corresponds to perturbations in an inhationary stage [3]. To get the factor 
A, we have to match the value of the short wavelength perturbations to that of the 
long wavelength. The result is: 


A 


H 


j.3/2 


( 2 . 21 ) 


with H being evaluated at rj rNj -^ Cs and k are the speed of sound and the wave 

Csk 

number, respectively. These results are in agreement with [3] 
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Chapter 3 
Conclusion 


In this paper, a Quintessential inflation scenario from the Mimetic Dark Matter 
model [T] has been presented. The potential used to produce such a scenario is 
defined on three time intervals, one during inflation(f = 10“^® — one after 

inflation that extends for and the third after that. The parameters of the 

potential were set in a way to produce 70 e-folds inflation and to have an energy 
density corresponding to the one measured today, representing Dark Energy. The 
scale factor after inflation is that of an accelerating universe, in contrast to a decel¬ 
erating De-Sitter Universe as it has been presented in the literature [TO] . 

However, in the non-modifled General theory of Relativity, an energy density of 
matter/radiation is accompanied by a decelerating Universe, this is why usually it 
is required to have a decelerating universe after inflation. But in this case, from the 
equations of MDM, one obtains an energy density of matter/radiation dominated 
universe. The important thing is the energy density rather than the scale factor. 
Since in this model we got the required energy density but not the ’’usual” scale 
factor, we can avoid the problem of explaining why the universe should decelerate 
after inflation. Rather, at the end of inflation, the universe looses enough energy 
for it to remain accelerating, but with a slower rate than the one during inflation. 
Moreover, concerning the structure of the potential, this form has been presented in 
physical problems other than Cosmology, mainly in electrostatics na, in order to 
produce physically acceptable and non-diverging electric fields. Any discontinuities 
or divergences in the potential should not be considered dangerous, as long as the 
Physical quantities are smooth and well behaved. Of course we still have to check 
whether this will still result in the required Nucleosynthesis and we have to see if 
the temperature perturbations that arise matches those of the CMB. These will be 
handled in future work. 
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